Invisible Decay of Muonium: 
Tests of the Standard Model and Searches for New Physics 
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The invisible decay of muonium atom M —> invisible, M — fi + e~ bound state, is expected 
to be a very rare process which has never been experimentally tested. In the standard model it 
goes via the n + e~ annihilation into two neutrinos, with the branching fraction predicted to be 
Br(M — } v^Ve) = 6.6 x 1CP 12 with respect to the ordinary muon decay rate. Using the reported 
experimental results on searches for muonium-antimuonium conversion at PSI we set the first limit 
Br(M — > invisible) < 4.6 x 10~ 3 (90%C.L.), while still leaving a big gap of about nine orders 
of magnitude between this bound and the predictions. To improve substantially the limit, we 
proposed to perform an experiment dedicated to the sensitive search for the M —¥ invisible decay. 
A feasibility study of the experimental setup shows that the sensitivity of the search for this decay 
mode in branching fraction Br(M — > invisible) at the level of 10 -12 could be achieved. If the 
proposed search results in a substantially higher branching fraction than predicted, say Br(M — > 



invisible) ~ 10 



this would unambiguously indicate the presence of new physics. We point out 



that such a possibility may occur due the muonium transition into a hidden sector and consider, as an 
example, muonium- mirror muonium conversion in the mirror matter model. A result in agreement 
with the Standard Model prediction would provide a theoretically clean check of the pure leptonic 
bound state annihilation through the charged current weak interactions, and provide constraints for 
further attempts beyond the Standard Model. We believe our work provides strong motivations to 
perform the proposed experiment on search for the invisible decay of muonium in near future. 



PACS numbers: 14.80.-j, 12.20.Fv, 13.20.Cz 



INTRODUCTION 



Experimental studies of particles invisible decays, i.e. 
transitions to an experimentally unobservable final state, 
played an important role both in development of the 
Standard Model (SM) and in constraining its extensions 
[l|. It is worth to remember the determination of the 
number of lepton families in the SM through the preci- 
sion measurements of the Z — > invisible decay rate. In 
recent years, searches for invisible particle decays have re- 
ceived considerable attention. One could mention exper- 
iments looking for extra dimensions with invisible decay 
of positronium (Ps=e + e~ bound state) [2], |3j|, baryonic 
number violation with nucleon disappearance at SNO Q , 
BOREXINO [|, and KamLAND [(| , see also Q, elec- 
tric charge nonconserving electron decays e~ — > invisible 
Q, neutron-mirror neutron oscillations at PSI Q and 
the ILL react or 10l| , neutron disappearance into another 
brane world [11| . and motivated by various models of 
physics beyond the SM searches for invisible decays of 
7T° mesons at E949 (H], V and '/mesons at BES [l|[, 
heavy B meson decays at Belle Q, BaBAR [Hj], BES 
[l6| and invisible decays of the Upsilon(lS) resonance at 
CLEO [17]. There are also proposals for new experiments 
to search for electric c harg e nonconservation in the muon 
decay u + —± invisible [181 ], and mirror type dark matter 
through the invisible decays of orthopositronium in vac- 
uum [191 ]. 

Muonium (M), the /i + e~-bound state, is a particu- 
larly interesting system for high precision tests of the 



SM and searches for new physics. Many interesting ex- 
periments performed or planned with muonium were mo- 
tivated by tests of bound state quantum electrodynamics 
(QED) in measurements of the muonium hyperfine split- 
ting [20| and ls-2s interval [21[, searches for the lepton 
number violation in muonium to anti-muonium conver- 
sion [13] , tests of fundamental symmetries, such as CPT 
[23| , probe of antimatter gravity in the free gravitational 
fall of muonium [24[ , and other areas of research [2514271 ] . 

As far as the muonium invisible decay is concerned, 
this process has never been experimentally tested. How- 
ever, there are several interesting motivations for such 
experiment to be performed. First, the decay is predicted 
to exist in the SM at the experimentally achievable today 
level. Hence, its observation for the first time would be an 
interesting test of the SM. Second, the decay may occur 
in some low-mass dark matter scenarios, most of which 
require coupling between the SM and hidden sectors. For 
instance, we show that in the mirror matter model such 
coupling could significantly enhance the M — > invisible 
decay rate, thus making it very attractive for direct high 
sensitivity searches. If the M — > invisible decay is ob- 
served at a rate higher than the SM prediction, it would 
be a strong evidence for the existence of new physics be- 
yond the SM. 

In this paper we obtain the first limit on the de- 
cay M — > invisible and show that it could be signifi- 
cantly improved in a new proposed high-statistics and 
low-background experiment. We also show that the ex- 
pected level of the experimental sensitivity allows for the 
observation of this decay mode for the first time at a rate 
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predicted by the SM. The rest of the paper is organized 
as follows. In Sec. II we briefly review the standard 
model decay rate of muonium and phenomenology of the 
muonium to anti-muonium conversion. In Sec. Ill the 
exact mirror model, the effect of oscillation of ordinary 
muonium to the mirror one, and its experimental conse- 
quences are discussed. The first limit of the branching 
fraction for the decay M —¥ invisible is obtained from 
available experimental data in Sec. IV. The experimental 
technique and the preliminary design of the experiment, 
detector components, simulations of the signal and back- 
ground sources, as well as the expected sensitivity are 
discussed in detail in Sec.V. Sec. VI contains concluding 
remarks. 



II. MUONIUM DECAY IN THE STANDARD 
MODEL AND PHENOMENOLOGY OF M -> M 
CONVERSION 



The muonium atom consits of a positive muon and 
an electron, which are leptons from two different genera- 
tions. To our current knowledge these are particles with- 
out any known internal structure. This makes muonium 
an ideal system for testing QED and fundamental sym- 
metries in physics [281 ] . and allows to calculate muonium 
properties to very high accuracy within the framework 
of the bound state QED. For example, for the hyperfine 
structure the theoretical predictions and measurements 
agree substantially better than for hydrogen atoms [29j . 

As discussed previuosly, muonium atom, similar to 
the lightest known exotic hydrogen-like atom positron- 
ium, is bounded by the electromagnetic (e-m) interac- 
tion. However, differently from positronium the muo- 
nium cannot self-annihilate through the e-m interaction 
because it would violate the lepton number conservation. 
Instead, the SM allows the self-annihilation of muonium 
into neutrino pair through the lepton number conserving 
weak interaction. At the current level of experimental 
and theoretical precision this are the only interactions 
present in the muonium system . 

In the SM muonium is unstable mainly due to the de- 
cay of muon fi + — > e + v e v li . Its decay rate in vacuum 
coincides with those of the muon decay given by 
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(1) 



Here Gf = 10 m~ is the Fermi constant and m p , 
are the proton and muon masses, respectively. In mat- 
ter, muonium is typically formed in the singlet or triplet 
state, with the total angular momentum equal or I, 
respectively. The SM predicts direct annihilation of the 
triplet muonium, J PC — 1 bound state, into neutrino 
antineutrino pair M — ¥ v e D^ with a very small decay 
rate. The corresponding branching fraction Br{M — ¥ 



Br(M — » invisible) — 
= 4W(- 



r(M -» v e v^) = 

^) 3 «6.6x 10- 12 , 



(2) 



where a is the fine-structure constant, and m e is the elec- 
tron mass. The singlet muonium cannot decay into two 
neutrinos, as it contradicts to momentum and angular 
momentum conservation simultaneously. 

Very interesting fearure of muonium is the possibil- 
ity of conversion(or oscillation) into its antiatom, i.e. 



the /i~e + bound state [311 - 4331 ] . This reaction violates 
the conservation of lepton flavor numbers by two units 
(AL e ,/i = ±2) that makes searches for the M — M con- 
version especially interesting due to experimental obser- 
vations of lepton flavor violation in neutrino oscillations. 

The existence of muonium to anti-muonium conversion 
is predicted in different extensions of the SM. The sim- 
plest way to understand the phenomenology of M — M 
conversion is the use of the effective four fermion inter- 
action of the (V — A)(V — A) type 32j, namely 
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where G MM is a coupling constant characterizing the 
strength of a new flavor violating interaction. In the 
absence of an external magnetic field the muonium and 
anti-muonium have the same ground-state energy levels. 
Flavor violating interaction (3) would cause a splitting of 
their energy levels by amount [32|, l33| | 
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Here n is the principal quantum number of the muonium 
atom and an = T " e+T " M is the Bohr radius of the muo- 
nium. Numerically, for the ground state of muonium 
(n = I) 



S = 1.5 x 10- i2 { 
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The M — M conversion is analogous to the K° — K° 
mixing. If a muonium atom is formed at t = in vacuum, 
it could oscillate into an anti-muonium atom. For a small 
t value the probability of the oscillation is represented in 
the form |32j ] 
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Here T M = T(fi + — > e + v e D^) is the muon decay width. 
The total conversion probability after integration over 
time is equal to 
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The best current experimental limit on the M — > M' 
conversion leads to bound \G m m\ < 0.003 • Gp [23 ■ It 
should be noted that in the presence of collisions with 
residual gas molecules or Zeeman splitting due to nonzero 
external magnetic field, the M — M transitions become 
highly suppressed. 
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III. MUONIUM CONVERSION IN MIRROR 
MATTER MODEL 

The idea that along with the ordinary matter may exist 
its exact mirror copy is an old one [3~il l35j . This new hid- 
den gauge sector is predicted to exist if parity is unbro- 
ken symmetry of nature, for an excellent review see Ref. 
36] . In accordance with this idea each ordinary particle 
of the SM has a corresponding mirror partner of exactly 
the same mass as the ordinary one, so that the properties 
of the mirror particles completely mirror those of the or- 
dinary particles. The parity symmetry interchanges the 
ordinary particles with the mirror ones. For example, the 
mirror proton and mirror electron are stable and inter- 
act with the mirror photon in the same way in which the 
ordinary proton and electron interacts with the ordinary 
photons. The mirror particles couple very weakly to the 
ordinary particles, hence they are not produced abun- 
dantly in our experiments. In the modern language of 
gauge theories, the mirror particles are all singlets under 
the standard G = SU(3) <g> SU{2) L ® U(1) Y gauge inter- 
actions [37H39( . Instead the mirror particles interact with 
a set of mirror gauge particles, so that the gauge symme- 
try of the theory is doubled, i.e. the ordinary particles 
are singlets under the mirror gauge symmetry. Parity is 
conserved because the mirror particles experience V+A, 
i.e. right-handed mirror weak interactions while the ordi- 
nary particles experience the usual V-A, i.e. left-handed 
weak interactions. Mirror matter, which interacts with 
ordinary matter presumably via gravity and possibly by 
other very weak forces, is dark in terms of the SM inter- 
actions, and could be a viable candidate for dark matter, 
see e.g. Ref s . [39M43| . For instance, it is argued that an- 
nual modulations of the signal observed by the DAMA 
Collaboration are caused by the mirror dark matter scat- 



M+ 



tering in their detector [4J-|47|. 

The gauge group of our world and mirror world is as- 
sumed to be 133 



Gtot = Gsm <& Gm, 



(8) 



where the SM gauge group Gsm coincides with the mir- 
ror world gauge group Gm- The interaction between our 
and mirror sectors could be transmitted by some gauge 
singlet particles interacting with both sectors. Such kind 
of interaction could explain the baryon asymmetry of the 
Universe [HI, some fraction of dark matter in the Uni- 
verse [43|], and can also induce the particle mixing and 
oscillation phenomena between the ordinary and mirror 
sectors. Any neutral, elementary or composite parti- 
cle, can have a mixing with its mirror duplicate, such 
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FIG. 1: The double degeneracy between mass eigenstates of 
ordinary (M) and mirror (M') muonium is broken when a 
small mixing is included. 



as photon-mirror photon [49i . l50J, neutrino- mirror neu- 
trino [33, HH, etc. which results in testable experimental 
effects. For example, the neutron-mirror neutron mix- 
ing via a small mass term e(nn' + n'n) proposed in [48| . 
results in ordinary neutron anomalous disappearance, in 
addition to their decays or absorption due to SM inter- 
actions. At present, there are performed and proposed 
searches for mirror matter via the invisible decay of or- 
thopositronium in vacuum through neutron- 

mirror neutron oscillations I48l. |56l l57j . and via Higgs- 
mirror Higgs mixing at LHC [58K60J. 

As mentioned previously, in mirror world model there 
must exist the mirror muonium which is the bound state 
of a mirror muon // and a mirror positron e' + with the 
same mass and decay width as the ordinary muonium. 
We could assume the existence of a (super)weak interac- 
tion invariant under the gauge group Gtot-, which allows 
transitions between the ordinary and mirror muonium. 
Such effective interaction can be written in a form anal- 
ogous to Q 
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where Gmm' is a coupling constant characterizing the 
strength of the M — M' transition, and e' and \i' are the 
mirror electron and muon fields, correspondingly. The 
interaction ^) leads to conversion of ordinary muonium 
to mirror muonium, as schematically illustrated in Fig. 
[1] The interaction © breaks the degeneracy between 
M and M' states, so that the vacuum energy eigenstates 
are M+ = (M + M')/\/2 and M— = (M - M')/V2, 
which arc split in energy by AE given by an expression 
similar to ((4]), see Fig[TJ The interaction eigenstates are 
maximal combinations of mass eigenstates which implies 
that M oscillates into M' . Thus, a system which is pure 
muonium at t = will develop an admixture of mirror 
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muonium at a later time. In closed analogy with the 
case of muonium anti-muonium conversion one can find 
that the probability of seeing the system in vacuum decay 



ordinary muonium M 



P(M') = - 



e fi 

e + v e v ll + e~ is given by 
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where A = 0.45 x 10 6 sec -1 , or 3 x 10~ 10 eV, and A is any 
additional splitting of M and M', by external electromag- 
netic fields. A detail discussions of the effects of collisions 
and external fields on oscillation probability of a similar 
system, positronium, can be found in Ref.|6l|. Estimat- 
ing 5 by using (9) results in the integral probability of 
muonium to mirror muonium conversion determined by 
(for A = 0) 
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Because mirror muonium decays into a mirror elec- 
tron, positron and neutrinos the experimental signa- 
ture of the M — M' conversion is the invisible decay 
M — > invisible of the ordinary muonium in vacuum. 
Current bounds on coupling constant Guw are rather 
weak. For Gmm 1 = Gp the branching fraction of the 
ordinary muonium decay into invisible mirror state is 
Br(M — > invisible) — 2.56 x 10~ 5 , which is seven or- 
ders of magnitude higher than the predicted by the SM 
branching fraction of Eq.(2). Thus, we see that in the 
mirror model it is possible to have nonzero conversion of 
our muonium to mirror muonium. The signature of such 
conversion is the decay M — ¥ invisible and, moreover, 
it is possible to expect some enhancement of this decay 
rate. Similar to the M — M conversion, the probability 
Pmm 1 can be affected by an additional splitting of M 
and M' states due to an external electric or magnetic 
field [H,[6l|. It might also be suppress, if there is a high 
collision rate of muonium atoms with the cavity walls or 
residual gas molecules in the experiment. 

Note that some extensions of the SM allow the M — > 
invisible decay. For instance, in model with additional 
sterile neutrino, interaction 

L = G'pj^l - j 5 )eD s j u (l - j 5 )i/ a + h.c. 

results in invisible muonium decays into sterile neutrino 
M — > v a v a . However, constraints obtained from the 
agreement between the measured and predicted proper- 
ties of the ^-decay [l[ lead to a strong bound T{M — > 
v s v a < O(10 _2 )r(M — > v e v^) on sterile neutrino decay 
width, which makes it not very exciting for further con- 
siderations. 



data. If an exotic M — > invisible decay exists, it would 
contribute to the total muonium decay rate: 



' M 



T M (M all) = r M +r(M ->• invisible)^... (12) 



and, hence decrease the determined muonium lifetime 
tm in vacuum. To estimate the allowed contribution 
of T(M — > invisible) to Eg . (jT2"j) . one could compare 
the measured muonium decay rate in vacuum to a pre- 
dicted muon decay rate, obtained in measurements which 
are not affected by the muonium disappearance effect 
discussed above. To avoid confusion with the decay 
/i + — > invisible, here we assume that both /x + and e~ 
disappear in the detector. 

In order to obtain limit on the branching fraction 

Br(M — > invisible) — T(M — > invisible) /T(p — > evv) 

(13) 

we use the results of the most precise experiment on 
M — M conversion performed at PSI (22|. In this ex- 
periment muonium was formed in a layer of Si02 powder 
and escaped into the vacuum after diffusion to the layer 
surface with a diffusion length L = ^Dmt^, where Dm 
and tm are the muonium diffusion constant and lifetime, 
respectively [62j . The distribution of time intervals be- 
tween the stop of a muon and the detected signal from 
the muonium decay in vacuum is shown in Fig. [5j The 
number NJf c (t) of M decays observed as a function of 
time and can be well approximated by [62] 



N V M C {t) = fexp(-t/T M )VDMt 



(14) 



where / is a normalization constant. The fit of the ex- 
perimental distribution shown in Fig. [5] to the function 
(O results in 



r^ p = 2.247 ±0.047^s. (15) 

The value of the positive muon mean lifetime is taken 
from results of its measurements by the MuLan collabo- 
ration has reported to aprecision of 11 ppm |63j |. Using 
the new world average [l| 

t m = 2.197019(21) fis (16) 

and comparing Eq. (fT5|l and Eq. (fT6|) . one finds 

Br(M — > invisible) = 

T(M invisible) 



r(^ 



ill) 



< 4.6 x 1Q- 6 (90%C.L.). (17) 



There is still nine orders of magnitudes difference be- 
tween the limit of Eq. dTTl) and the SM prediction (2). In 
the next Section we show how this limit can be signifi- 
cantly improved in the new proposed experiment. 



IV. EXPERIMENTAL LIMIT ON THE 

M -> invisible DECAY 

Consider now bound on the invisible decay of the M 
state, which can be obtained from existing experimental 



V. DIRECT EXPERIMENTAL SEARCH FOR 

THE M -> invisible DECAY 

The decays M — > invisible are rare events and their 
observation presents a challenge for the detector design 
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FIG. 2: The distribution of time intervals t between the stop 
of a muon in the SiC>2 target and signal from the decay of 
muonium escaping into vacuum for t > 2000 ns in the PSI 
experiment on search for M — M conversion [22] ]. Experi- 
mental points and dashed curve representing an exponential 
background are from Ref.[23|. the solid line is our fit to the 
data. 



and performance. Here, we focus mainly on discus- 
sions of the experimental setup to search for the decay 
M — > invisible in vacuum, which is also sensitive to 
the muonium - mirror muonium conversion. The similar 
setup without vacuum requirements is simpler, it would 
provide better sensitivity and can be used for the first 
observation of the decay M — > invisible and search for 
new exotic channels of this decay mode which are not 
affected by the presence of matter of external fields. 

The main components of the experimental setup to 
search for the invisible decay of muonium are schemat- 
ically illustrated in Fig. O see also [18j. The setup is 
equipped with a high efficiency muon tagging system, 
high hermeticity electromagnetic calorimeter (ECAL), 
and an intelligent trigger system. The experiment em- 
ploys a surface /i + beam, which is produced in a target 
and transported to the detector in an evacuated beam 
line tuned to ~ 26 MeV/c. Such the world's brightest 
continuous surface muon beam with intensity ~ 10 7 /i/s 
is available at the Paul Scherrer Institut (PSI) 64] . This 
beam was used, for example for a sensitive search for 
M — M conversion [22|. Positively charged muons pass 
through a ~ 100 /im thick beam counters (Si. 2) and fo- 
cused into a vacuum cavity through a narrow aperture 
closed by the beam counter S3 and after passing through 
the counter S4 strike the SiC>2 aerogel (or SiC>2 pow- 
der) target (T) used for the muonium atom formation 
[60I |66| . The energy of entering muons is degraded by 
the counters material to maximize the muon stopping 



rate. Muonium atoms are formed by the electron cap- 
ture with efficiency ~60 % per fi + stopped in the target. 
Most of the atoms emerge from the target grains into 
the intergranular voids. With a mean-free path of ~ 1 
cm, muonium is able to diffuse through the network of 
voids over distances longer than the target thickness and 
escape through the surface into vacuum [28j . Muoniums 
undergo of the order 10 5 — 10 7 collisions with the silica 
grain walls, those number depends on on the depth of 
muonium formation, and approach thermal equilibrium. 
Then, on average, 3.3% of them leave the target surface 
with thermal Maxwell-Boltzmann velocity distribution at 
the temperature of the target [65[. The fact that muo- 
nium confined to the voids is expected to be almost fully 
thermalized, was confirmed by a separate experiment on 
M's emitted into vacuum from a mesaporous silica film at 
cryogenic temperatures [67| ■ Alghough the M kinetic en- 
ergy distribution is nominally that of Maxwell-Boltzman 
emission, one might expect a higher-energy tail of M's 
formed from backscattered muons that never approached 
thermalization. These M events are distinct from the 
thermal M's that has diffused out of the target, however 
their intensity is expected to be very small. The fraction 
of muonium atoms produced in the target that decay ei- 
ther in the target or in vacuum can be determined relative 
to the number of muons on the target with a technique 
described in Ref . [22| . 

The target is surrounded by a hermetic 4ir electro- 
magnetic calorimeter (ECAL) to detect energy deposi- 
tion from the decay M — > all of muoniums produced 
in T. As shown in Fig. [31 before muons reach the en- 
trance to the vacuum cavity, they bend in a magnetic 
field. The purpose of employing the magnet is twofold: 
(i) to provide a transverse kick to positive muons in or- 
der allow them to enter the vacuum cavity through the 
narrow aperture, and (ii) to detect photons, positrons, 
or muons that could escape cavity through the entrance 
aperture by a set of ECAL counters placed around the 
muon bend region. This additional detector is placed 
upstream of the entrance aperture, as shown in Fig. [3j 

The energy deposition readout in the ECAL is trig- 
gered by a tag signal of the muon appearance on the 
target, which is defined as the coincidence of the four 
signals from a muon passing the beam counters Si— 4. To 
enhance significance of the muon tag the time-of-flight 
information can be used. For the muon beam momen- 
tum of 26 MeV/c , the latter corresponds to about 1.3 
ns per 10 cm of the muon path length. 

For the ordinary muonium decay ([T]) the experimen- 
tal signature is the ECAL energy deposition from a sin- 
gle decay positron with energy E e + = — E Ve — , 
where E Vc , E v are the electron and muon neutrino en- 
ergy, respectively. The experimental signature of the 
M — > invisible decay is the apparent disappearance of 
the energy deposition E = + m e in the ECAL. In 
other words, the signature of the M — > invisible decays 
is an event with the sum of the ECAL crystal energies 
deposited by the final-state particles equal zero. Zero en- 
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FIG. 3: Schematic illustration of the experimental setup to 
search for the muonium invisible decay. The beam of surface 
fi +, s passing through the beam defining counters S1-4 is fo- 
cused into a vacuum cavity through a narrow aperture and 
strike the Si02 aerogel target (T) used for the muonium atom 
formation. Shown are also the 4n hermetic BGO electromag- 
netic calorimeter (ECAL), the ECAL endcap counter EC used 
as a light guide for the light produced in the beam counter 
S4, and the magnet used to deflect the beam. The counters 
S1-3 and the upstream ECAL counters are also used as a veto 
against photons, decay positrons or backscattered muons that 
could escape the cavity through the entrance aperture. 



ergy is defined in this case as an energy deposition below 
of a certain ECAL energy threshold, E to t < E t h- The 
expected distribution of energy deposited in the ECAL 
from /i + 's stopped in the target is shown in Fig. fJJ The 
distribution is a sum of two spectra from /i + — » all and 
M —¥ all decays and is discussed in detail below in Sec. A. 

To estimate the sensitivity of the proposed experiment 
a feasibility study based on GEANT4 [n| Monte Carlo 
simulations have been performed. The beam of positive 
muons is stopped in the central part of the cylindrical 
target of S1O2 aerogel with the density about 35 mg/cm 3 
and of thickness 8 mg/cm 2 and supported in vacuum by 
an aluminum foil with an inclination with respect to the 
muon beam axis. Similar target was previously used as a 
convertor of muon to muonium atoms in the experiment 
of Ref.[H|. The ECAL is an array of ~ 100 BGO coun- 
ters each of 52 mm in diameter and 220 mm long, which 
was previously used in the PSI experiment on precise 
measurements of the n — ¥ e + v decay rate [69| . Timing 
and energy deposition information from each BGO crys- 
tal can be digitized for each event. The processingof the 
BGO counter signals is described in detail in Ref . |69| . see 
also 0, 0| . 



A. Background for the M —¥ invisible decay 

The background processes for the M — ¥ invisible 
decay can be classified as being due to beam-related, 
physical, and detector-related backgrounds. To inves- 
tigate these backgrounds down to the level Br(u —¥ 
invisible) < 10 -12 with the full detector simulation 



LU 

< 




10 



10 



10 



MeV 



FIG. 4: The expected distribution of energy deposition in 
the ECAL from 8 x 10 12 muons stopped in the target, cor- 
responding to the decays — ¥ all and M — ¥ all. The 
peak around 1 MeV corresponds to energy deposition from 
the e + e~ — > 27,37 annihilation of decay positrons stopped 
in the vacuum cavity. The arrow shows the energy threshold 
for the decay M —¥ invisible detection. The dashed curve 
represents the signal from the decay M —¥ invisible if it exist 
at the level predicted by the SM. 



would require the generation of a very large number of 
muon decays resulting in a prohibitively large amount of 
computer time. Consequently, only the most dangerous 
background processes are considered and estimated with 
a smaller statistics combined with numerical calculations. 

The beam-related backgrounds produce the fake muon 
tag and can be categorized as being due to a beam par- 
ticle misidentified as a muon, or several beam particles 
which produce fake muon tag due to accidental coinci- 
dence of signals from Si_4. The first type of background 
occurs, e.g. due to the production of slow protons in 
the target, which enter the detector and produce zero 
decay energy. Incoming neutrons could scatter in the 
Si_4 and being accidentally misidentified as /i + could 
also contributed to the beam background. Identification 
of the incoming particle as a muon based on the require- 
ments of the delayed by the muon time-of-flight coinci- 
dence between the beam counter signals suppressed the 
single-beam background down to the level < 10 -13 . This 
estimate is obtained under assumption of having Gaus- 
sian shape with the time resolution of ~ 1 ns for the 
distributions of time-of-flight between counters Si -4. It 
is also assumed that the admixture of the other charged 
particles in the beam is below 1%, which although de- 
pends on a particular experimental environment. 

For the design shown in Fig|31 the required efficiency 
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FIG. 5: The fraction of events from the decay M — > e + v^u e + 
e~ with the positron energy E < Eth as a function of Eth- 



for the M decay energy detection can be obtained only 
by keeping the amount of passive material in the region 
of vacuum cavity as small as possible. For example, to 
remove dead materials from the vacuum cavity walls the 
cavity could be made directly in a big single crystal or 
out of a few ECAL central crystals. The light signals 
produced in the S4 scintillator counter could be readout 
through the SiC>2 transparent target and the ECAL end- 
cap crystal (EC) which acts as a light guide, see Fig. 
[H The S4 signals could be distinguished from the EC 
signals due to their significantly different decay times by 
using the technique described in detail in Ref. 3]. In the 
presented simulations we did not consider too compli- 
cated design of the setup and try to keep it as realis- 
tic as possible. The reported further analysis takes into 
account active materials of the ECAL, passive materials 
from the target, vacuum cavity walls, and from the ECAL 
crystals and the target wrapping. The following main 
sources of physical- and detector-related backgrounds are 
identified and evaluated: 

• the principal muon decay y, — > evv into the final 
state with the positron kinetic energy Eki n less 
than the detection energy threshold By, ( ~ 100 
keV). Indeed, if Ekin < E t h the event becomes in- 
visible. To suppress this background, one has to use 
as low as possible threshold E t h and to performed 
the experiment with a well separated positive muon 
beam with an extremely small contamination of 
negative pions or muons which could mimic the true 
signal. However, even if the positive muon decays 
into a low energy positron that stop in the cavity, 
the later would annihilate into two (or three) pho- 
tons at a lifetime scale of the order of a few ns. 



Thus, for such events the minimum energy depo- 
sition in the ECAL will be m e + + m e - ~ 1 MeV, 
i.e. well above the threshold, making these events 
visible, see Fig. 2) 

Another way to lose the decay energy is due to 
the annihilation gammas photo- absorption and/or 
Compton scattering in the target. In this case, 
when almost all annihilation energy is deposited 
in T the event becomes invisible, which results in a 
fake M — > invisible signal. To suppress this back- 
ground the target should be optimized in size and 
made of a low-Z material to minimize cross-section 
of the photo-absorption which is a p h a ~ Z 5 . For 
example, for a target made of plastic scintillator 
the probability of both 511 keV photons energy ab- 
sorption in a volume of ~ 1 cm 3 is found to be less 
than 10 -8 [18j. In the Si02 target with the den- 
sity 35 mg/cm 3 the effect is smaller. The ECAL 
efficency with respect to detection of energy from 
the positron annihilation was checked in the exper- 
iment [3J on search for the Ps — > invisible decay. 
For events corresponding to 2, 37 annihilation of 
e + e~ pairs at rest with the ECAL energy depo- 
sition ~ 1 MeV, the upper limit on the branch- 
ing fraction of the reaction e + e~ — ¥ invisible was 
found to be Br(e + e~ — > invisible) < 10 -8 at 90% 
C.L. for the ECAL energy threshold of 80 keV. 

The absorption of annihilation photons in the cav- 
ity materials has been studied in the proposal on 
search for the oPs — > invisible decay in vacuum 
of Ref. [IH . Simulations show that the main con- 
tribution to the 7 inefficiency comes from the to- 
tal (due to photo-absorption) or fractional (due to 
Compton effect) photon energy loss in the material 
of the vacuum cavity. To suppress this background 
the cavity should be made of a low-Z material 
to minimize cross-section of the photo-absorption. 
Distributions of the energy deposited in the dead 
material surrounding the target region from annihi- 
lation events in the target were obtained with sim- 
ulations for a 0.84 mm thick aluminum pipe and a 
composition pipe made of 0.04 mm aluminum and 
0.800 mm carbon. For the later case the fraction 
of simulated 27 events with the energy absorbed in 
the cavity walls > 900 keV, i.e. energy deposited in 
the ECAL is E < 100 keV, was found to be < 10~ 8 . 

In Fig. [5J the partial muon decay rate AT^ into a 
positron with E e + < Eth is shown as a function of 
E t h- Taking into account that energy of positrons 
that stop in the cavity is typically Ekin < 2 — 3 
MeV, the fraction of such e + 's is estimated to be 
Pr < 10~ 4 . Combined probability to get energy 
deposition in the ECAL from positron annihilation 
in the cavity < 100 keV is estimated to be P27 ^ 2 x 
10~ 8 . Therefore, this background from inefficient 
detection of low energy positrons allows potentially 
to reach sensitivity in the branching ratio of the 
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invisible muonium decay as small as Pr ■ -P27 — 
2 x 10~ 12 , assuming the detection energy threshold 
is as low as E t h — 100 keV. 

• the loss of the muon decay energy in rare processes 
of energetic positrons e + + A — > invisible with an 
invisible final state could be induced either by elec- 
tromagnetic or weak interactions of the positron. 
For example, e + could lose almost all its energy due 
to emission of a hard photon in the bremsstrahlung 
process in the target or in the ECAL. The pho- 
ton could either penetrate the calorimeter with- 
out interactions, or could be photo-absorbed by an 
atomic nucleus resulting in the invisible final state 
consisted of secondary neutrons. However, due to 
the charge conservation, there is always a low en- 
ergy positron in the final state, which produce ~ 
1 MeV energy through the e + e~ annihilation of 
the positron at rest, thus making the event visi- 
ble. Combined analysis results in this background 
level < 10 -13 . The background from an energetic 
positron conversion into proton through the reac- 
tion e + + n — > p + v e induced by the charge current 
weak interaction is found to be negligible. 

• Another possible background could be due to the 
excitation of a long-lived nuclear state via the ra- 
diationless annihilation of an energetic positron 
with a K-shell electron. This is a 3-body reaction 
e + + e~ + A — > A*, where the e + e~ annihilation 
energy is absorbed by the nucleus A. The cross 
section for such reaction has not yet been stud- 
ied in details for the wide class of nuclear isotopes 
and full range of positron energies. By using the 
available upper limit on the resonant cross section 
a e + < 4.3 x 1(T 26 cm 2 at 99% C.L. obtained for 
isotope 115 In with a mono-energetic positron beam 
of about 90 keV kinetic energy [7l| we estimate this 
background to be < 10 -13 , assuming the 115 In con- 
tamination in the cavity and target materials to be 
at the level below 1 ppm. More detailed study of 
this background source is required. Note, that in 
principle, it is possible to excite a nucleus long-lived 
state with a lifetime r > 60/is. However, such ex- 
citation levels are present in specific isotopes, such 
as 115 In, those admixture is expected to be small. 

• incomplete ECAL hermiticity: our study identified 
possible background to the signal as due to ener- 
getic decay positrons escaping the detection region 
though the cavity entrance aperture. This effect 
increase the disappearance rate of muonium and 
therefore must be addressed. Consider, e.g. the 
case when a muon decays either in flight or in the 
target into a fast positron with momentum pointing 
exactly to the entrance aperture. Then, the decay 
positron could be undetected in beam counters S3 
and S2 due to their inefficiency. The same effect 
could occur if the incoming muon backscatters ei- 
ther in S4 or in the target without losing too much 



energy, and escapes the detection in counters 82,3. 
However, due to presence of the magnetic field in 
the vicinity of the entrance to the cavity the trajec- 
tory of the escaping positron or muon is bended up 
and it would be detected by the ECAL counters. 

The probability for a particle to escape detection 
in this case can be estimated as 

P esc ~ P a ■ P m ■ C2 • Ca • Cecal (18) 

where P a ,(2, C3i an d Cecal are, respectively, the 
probability for particle to pass through the entrance 
aperture, the inefficiencies of beam counters S2, S3 
and ECAL counters to detect the particle. 

To suppress this type of background the entrance 
aperture should be reduce to as much as possible 
size and should be closed by as high as possible ef- 
ficiency counters 82,3, as shown in Fig. [3J which 
act as the beam defining and also as the veto coun- 
ters. Then, the background could be suppressed by 
requiring an absence of activity in the beam coun- 
ters after detection of incoming muon. Assuming 
isotropic distribution of backscattered muons or de- 
cay positrons, inefficiency for particle detection in 
82,3 of ~ 10~ 2 — 10~ 3 , the diameter of the entrance 
aperture of ~1 cm, and inefficiency of the upstream 
ECAL detection ~ 10 -4 leads to the final suppres- 
sion of this source of background down to the level 
of at least ~ 10~ 13 . 

• the leak of muonium atoms through the entrance 
aperture into the region of lower detection effi- 
ciency could also contribute be the disappearance 
rate of muonium. However, assuming that muon- 
inums leaving the target are thermalized and have 
kinetic energy below eV (300 K), the effect is sup- 
pressed to a negligible level by closing the aperture 
with the counter S3, as shown in Fig. [3J 

In Table I contributions from the previously discussed 
background processes are summarized. The dominant 
background source is due to the absorption by pas- 
sive materials of photons from the annihilation of slow 
positrons in the cavity. To cross-check this background, 
we estimate its level in the signal region by using avail- 
able results from measurements of Ref.Q and the pro- 
posal on the search for oPs — > invisible decay in vacuum 
of Ref. [lH in a different way. In Fig. [4] the expected 
distribution of energy deposition in the BGO calorimeter 
from the decays of 8 x 10 12 /i + 's stopped in the target 
is shown. The spectrum represents the sum of /i + — > all 
and M — > all distributions . The part of the spectrum 
above > 1 MeV is calculated from the Michel spectrum 
convoluted with the ECAL (Gaussian) energy resolution. 
The peak around ~ 1 MeV is from the fraction of de- 
cay positrons (~ 10~ 4 ) with energy below of a few MeV 
that are stopped in the cavity, i.e. either in the target 
or in the cavity walls, and annihilate into 2 or 3 pho- 
tons. The low energy tail below 1 MeV is described by a 
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TABLE I: Espected contributions to the total level of back- 
ground from different background sources ( see text for de- 
tails). 



Source of background 


Expected level 


fake muon tag 


£ 10" 


■ra 


inefficiency of of slow positrons detection " 


< 2 x 10" 


■12 


e + + A — » invisible 


< 10" 


-13 


ECAL hermiticity 


~ 10" 


■13 


Total ( conservatively) 


~ 2.3 x 10" 


VI 



a The threshold for energy deposited in the ECAL from the decay 
e+'s annihilation is 100 keV. 



function f(E E +) — fi(E e +) + f2{E e +), which is a sum of 
two distributions of the annihilation energy in the ECAL 
normalized to the same number of positrons annihilated 
in the cavity. The function fi(E e +) is an experimen- 
tally measured distribution taken form the experiment 
on Ps — t invisible [3| for positrons annihilated in the 
Si02 target, which did not take into account the annihila- 
tion photon absorption in the cavity walls. The function 
f2 (E e + ) is taken from the proposal and corresponds 
to the simulated energy deposition in the ECAL minus 
energy absorbed in the cavity walls. The sum function 
f(E e +) is then extrapolated to zero energy resulting in 
prediction of about 8 ± 2 background events in the signal 
region for 8 x 10 12 /i + 's stop in the target, which is some- 
what smaller, but still in a reasonable agreement with the 
conservative number of about 18 events obtained from 
from Table 1. The error of the above estimate is denned 
by the uncertainty in the extrapolation procedure. 

B. Sensitivity of the proposed experiment 

The significance of the M — > invisible decay discovery 
with such detector, scales as [z^, [73| 

S = 2 ■ (Vn s + n b - v/n£) (19) 

with 

n s = n t± efBr(M —> invisible)t (20) 

and the branching ratio of the muonium invisible decay 
defined by 

Br(M -» invisible) = — (21) 

n^eft 

where and n s is the number of observed signal events 
(or the upper limit of the observed number of events), 
ni, is the number of background events, n M is the muon 
beam intensity, t is the experiment running time, e is the 
efficiency of the muonium production per incident muon, 
and factor / corresponds either to the total number of 
decayed M atoms (/ — 1), or to the fraction of M atoms 
that decay presumably in vacuum, not in the target (/ ~ 
0.033). 



Before defining the expected sensitivity, let us first dis- 
cuss several additional limitation factors. The first one is 
related to the relatively long muon lifetime and the cor- 
responding ECAL signal integration time. Indeed, to get 
the branching ratio Br{M — > invisible) ~ 10~ u , the 
ECAL gate duration r ff , and hence the dead-time per 
trigger, has to be 

T g > -t m x ln(Br(M -» invisible)) ~ 60 (is (22) 

in order to avoid background from the muon decays out- 
side the gate. The best sensitivity is expected at integra- 
tion gate t ~ 60 /is, however further more complicated 
analysis compromising the level of this background and 
increase of the pile-up noise might be necessary. The 
pile-up energy, which corresponds to energy deposited in 
the BGO ECAL by an additional undetected and uncor- 
related particle, increases values of the ECAL pedestals. 
The amount of additional energy in each BGO counter 
can be measured with the random trigger Q. In the Ps 
experiment [3j , for orthopositronium lifetime in the Si02 
aerogel target of 132 ns the ECAL gate duration rp s was 
chosen to be ~ 2 /is. This resulted in distribution of the 
sum of pedestals of all ECAL counters corresponded to 
the efficiency of "zero" signal detection as a function of 
the energy threshold. In order to keep the energy thresh- 
old as low as possible an algorithm to sum up the energy 
of all the ECAL crystals can be employed by exploiting 
the granularity of the calorimeter and fixing a zero en- 
ergy threshold for each individual crystal. Taking into 
account the ECAL granularity, the effective ECAL en- 
ergy threshold can be significantly reduced from 80 keV, 
used to define the signal range for the o — Ps — > invisible 
decay [1], to about 20 keV having the overall signal effi- 
ciency above 95% [70| . In the proposed experiment the 
longer gate will lead to an increase of the pile-up and pick- 
up electronic noise and hence to the overall broadening of 
the signal range, approximately by a factor \/t 9 /tp s ~ 5 
and, hence to an increase of the effective energy threshold 
roughly up to E th ~ 20 keVx5 ~ 100 keV. 

Another limitation factor is related to the dead time 
of Eq. (|2"2")l and, hence to the maximally allowed muon 
counting rate, which according to Eq. (|2"2"|) has to be 
< l/r g ~ 10 4 n + /s to avoid significant pile-up effect. To 
minimize dead time, the trigger electronics has to reduce 
the rate of events accompanied by an energy deposition 
in the ECAL. For example, one could implement a fast 
first-level trigger rejecting events with the ECAL energy 
deposition greater than E t h and, hence run the experi- 
ment at the rate ~ l/r M ~ 5 x 10 5 u/s. Assuming this 
rate, we anticipate 8 x 10 12 /i + on target and production 
of about 6 x 10 12 muonium atoms during 6 months of 
running time for the experiment. Out of them, about 
5.8 x 10 12 M's decay in the target, while about 2 x 10 11 
M's leave the target surface and decay in vacuum. For 
counting signal rate of ~ 10~ n per incident muon. as- 
suming beam intensity of ~ 5 x 10 5 u + /s at ~ 90% effi- 
ciency, it would require 1 week to accumulate one signal 
event. 
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In the background free experiment one could expect a 
sensitivity in the M — > invisible decay branching ratio 
of the order of 

Br(M -» invisible) < 10~ 12 , (23) 

assuming that in Eq. ([2~Tj) n s — 2.3. For M's that decay 
in vacuum, the sensitivity is 

Br vac (M -> invisible) < 10 -11 . (24) 

In the presence of background and in accordance with 
the SM prediction, the expected number of observed 
events in the signal region E < 100 keV is 

N M ~ 50 ± 7 events (25) 

out of which 18.4 events represent conservatively esti- 
mated overall background from Table 1. Taking into 
account (p~9|) . one can see that the observation of the 
M — > invisible decay with about 5 a significance could 
be possible. 

The statistical limit on the sensitivity of the proposed 
experiment to search for the decay M — > invisible due to 
transition into hidden sector is proportional to G 2 MM , and 
is set by its value, see ([TTj). Thus, to improve the sensi- 
tivity of (|24[) larger amount of muonium atoms decaying 
in vacuum is required. Therefore, the improvement of 
the efficiency for thermal M's production is crucial for 
further searches. 

Note, that in the case of the signal observation, to 
cross check the result, one could replace the target with 
another one of the same density, but not capable of muo- 
nium producing, and run the experiment with suppressed 
M decays, see e.g. (74[- I n this case the distribution 
of the energy deposition in the ECAL, shown in Fig|4] 
would contain mainly events from the decays fi + — > all 
and the signal from the decays M — > invisible should 
disappear. In the case of observation of a higher than 
predicted M — > invisible decay rate, there is another im- 
portant cross-check. Namely, as discussed in Sec. 3, one 
could slightly modify the experimental conditions with- 
out affecting the background , e.g. by increasing either 
the magnetic field in the cavity or the number of muo- 
nium collisions with residual gas molecules by increas- 
ing the gas pressure [H, |6l| . These would suppress the 
muonium-mirror muonium oscillations, and the observed 
signal should vanish. 

The performed analysis give an illustrative correct 
order of magnitude for the sensitivity of the proposed 
experiment. The simulations are performed without 
taking into account such effects as, e.g. pileup, and may 
be strengthened by more accurate and detailed Monte 
Carlo simulations of the concrete experimental setup. 



VI. CONCLUSION 

Due to its specific properties, muonium is an impor- 
tant and interesting probe of the SM and physics be- 



yond the SM both from the theoretical and experimental 
view points. In the SM the invisible decay of muonium 
atoms fi + e~ — > invisible into two neutrinos is expected 
to be a very rare process with the branching fraction 
predicted to be Br(M — > v^v e ) = 6.6 x 10~ 12 with re- 
spect to the ordinary muon decay rate. This process 
has never been experimentally tested. Using the re- 
ported experimental results on searches for muonium- 
antimuonium conversion at PSI we set the first limit 
Br(M — > invisible) < 4.6 x 10 -3 , while still leaving a 
big gap of about nine orders of magnitude between this 
bound and the predictions. 

To improve substantially the sensitivity, we proposed 
to perform an experiment dedicated to the search for 
the M — > invisible decay. The key point for the ex- 
periment is the presence of energy release from the an- 
nihilation of the low energy decay positrons in the de- 
tector. A feasibility study of the experimental setup 
show that the sensitivity of the search for this decay 
mode in branching fraction Br(M — > invisible) at the 
level of 10 12 could be achieved. Thus, the SM predic- 
tion for the M — > invisible decay to exist at the level 
of Br(M — > invisible) ~ 6.6 x 10~ 12 , could be experi- 
mentally tested for the first time. We point out that the 
M — > invisible decay rate could be enhanced by non- 
SM contributions. For instance, in the framework of the 
mirror matter model if the coupling strength between 
M and M' is large enough, say Gmm< ^ 10 _4 Gf, the 
decay M — > invisible could occur at a rate as high as 
the SM one. If the proposed search results in a substan- 
tially higher branching fraction than the SM predictions, 
say Br(M — > invisible) ~ 10~ 10 , this would unambigu- 
ously indicate the presence of new physics. A result in 
agreement with the SM prediction would provide a the- 
oretically clean check of the pure leptonic bound state 
annihilation through charged current weak interactions, 
and provide constraints for further attempts beyond the 
SM. 

The preliminary analysis shows that the quoted sen- 
sitivity could be obtained with a detector optimized for 
several its properties. Namely, i) the primary beam and 
the entrance aperture size, ii) the efficiency of the muo- 
nium production in the target and in vacuum, iii) the ma- 
terial composition and dimensions of the target and vac- 
uum cavity, iv) the efficiency of the veto counters Si -4, 
and v) the pile-up effect and zero-energy threshold in the 
ECAL are of importance. 

We believe our proposal, when paired with an existing 
BGO calorimeter, provides interesting motivations for 
the experiment on search for the M — > invisible decay 
to be performed in near future. This low-energy experi- 
ment might be a sensitive probe of new physics that is 
complementary to collider experiments. For example, it 
could also significantly improve the recently obtained 
modest bounds on the [i + — > invisible decay [II], push- 
ing it down to the region Br(u + — > invisible) ~ 10 -12 . 
A bound in this region will be of interest for several 
extensions of the Standard Model, see e.g. 75]. The 
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required high numbers of muonium atoms can be 
presently produced at PSI 22], or could be available 
from high intensity muon beams at future facilities such 
as the PRISM source at J- PARC ffjL the Project X at 
FNAL [73], or the neutrino factory |78j . 
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